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We study the scattering of massless scalar waves by a Kerr black hole by letting plane monochro- 
matic waves impinge on the black hole. We calculate the relevant scattering phase-shifts using the 
Priifer phase-function method, which is computationally efficient and reliable also for high frequen- 
cies and/or large values of the angular multipole indices (Z,m). We use the obtained phase-shifts 
and the partial-wave approach to determine differential cross sections and deflection functions. Re- 
sults for off-axis scattering (waves incident along directions misaligned with the black hole's rotation 
axis) are obtained for the first time. Inspection of the off-axis deflection functions reveals the same 
scattering phenomena as in Schwarzschild scattering. In particular, the cross sections are domi- 
nated by the glory effect and the forward (Coulomb) divergence due to the long-range nature of the 
gravitational fleld. In the rotating case the overall diffraction pattern is "frame-dragged" and as a 
result the glory maximum is not observed in the exact backward direction. We discuss the physical 
reason for this behaviour, and explain it in terms of the distinction between prograde and retrograde 
' motion in the Kerr gravitational field. Finally, we also discuss the possible influence of the so-called 

' superradiance effect on the scattered waves. 
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CN ; I. INTRODUCTION 

^ ■ . . 

^ . Diffraction of scattered waves provides the explanation for many of Nature's most beautiful phenomena, such as 
■ rainbows and glories. It has long been recognized that these optical phenomena have analogies in many other branches 
, of physics. They are of particular relevance to quantum physics, where plane wave "beams" are routinely used to 
probe the details of atoms, nuclei or molecules. Such experiments provide a deep understanding of the scatterer's 
\ physics and can be used as a powerful test of various theoretical models. The analogy can be extended also to 
gravitational physics and extreme astrophysical objects like black holes. In fact, black hole scattering has been the 
subject of a considerable amount of work carried out over the last 30 years (see for an extensive review). In the 
. case of astrophysical black holes it is unlikely that the various diffraction effects will ever be observed (although it is 
not entirely implausible that advances of current technology will eventually enable us to study interference effects in 
I , gravitationally lensed waves). However, it is nevertheless useful to have a detailed theoretical understanding of the 
scattering of waves from black holes. After all, a study of these problems provides a deeper insight into the physics 
of black holes as well as wave-propagation in curved spacetimes. 

The benchmark problem for black-hole scattering is massless scalar waves impinging on a Schwarzshild black hole. 
This problem is well understood [|[j8|, and it is known that it provides a beautiful example of the glory effect. Handler 
and Matzner ||^ have shown that the situation remains almost unchanged if, instead of scalar waves, one decides to 
cd ' "shoot" plane electromagnetic or gravitational waves towards the black hole. These authors have also considered 
on-axis scattering of gravitational waves in the case when the black hole is rotating Their results suggest that 
the scattering cross sections consist of essentially the same features as in the non-rotating case (there is a forward 
divergence due to the long-range nature of the gravitational field and a backward glory). In addition, they find some 
peculiar features that are, at the present time, not well understood. An explanation of these effects is complicated 
by the fact that they could be caused by several effects, the most important being the coupling between the black 
hole's spin and the spin/polarisation of the incident wave. Given that the available investigations have not been able 
to distinguish between these various effects, we feel that our current understanding is somewhat unsatisfactory. This 
feeling is enhanced by the fact that no results for the most realistic case, corresponding to off-axis incidence, have yet 
been obtained. 

This paper provides an attempt to further our understanding of the scattering from rotating black holes. Our aim is 
to isolate those scattering effects that are due to the spin of the black hole. In order to do this, we focus our attention 
on the scattering of massless scalar waves. For this case, the infalling waves have neither spin nor polarisation and 
therefore one would expect the scattered wave to have a simpler character than in the physically more relevant case 
of gravitational waves. However, one can be quite certain that the features discussed in this paper will be present also 
in the case of gravitational waves. It is, after all, well known that the propagation of various fields in a given black 
hole geometry is described by very similar wave equations. 
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Although we will re-examine the case of axially incident waves, our main attention will be on the more interesting 
off-axis scattering cross sections. These cross sections turn out to be quite different from the ones available in the 
literature. Obviously, they have two degrees of freedom (corresponding to the two angles and tp in Boyer-Lindquist 
coordinates). In addition we will show that the cross sections are asymmetric with respect to the incidence direction. 
In particular, the glory moves away from the backward direction as a result of rotational frame dragging that provides 
a distinction between prograde and retrograde motion in the Kerr geometry. 

We construct our differential cross sections using the well-known partial wave decomposition (for an introduction 



see |10 ) — the standard approach in quantum scattering theory. That this method is equally useful in black-hole 
scattering is well established We should point oui^, however, that alternatives (such as the complex-angular 
momentum approach (l^Jl^ and path- integral methods Jzl jl^Jl^ ) have also been succesfully applied to the black-hole 
case. In the partial wave picture all scattering information is contained in the radial wavefunction's phase shifts. 
The calculation of these phase-shifts must, apart from in exceptional cases like Coulomb scattering, be performed 
numerically. Various techniques have been developed for this task. Basically, one must be able to determine the phase- 
shifts accurately up to sufficiently large I partial waves that no interference effects are lost. This boils down to a need 
for many more multipoles to be studied as the frequency of the infalling wave is increased. In black-hole scattering 
several methods have been employed for the phase-shift calculation: Matzner and Ryan |^ numerically integrated 
the relevant radial wave equation (Teukolsky's equation). Since the desired solution is an oscillating function, this 
calculation becomes increasingly difficult (and time consuming) as the frequency is increased. Consequently, Matzner 
and Ryan restricted their study of electromagnetic and gravitational wave scattering to loM < 0.75 and I < 10. 
In order to avoid this difficulty. Handler and Matzner combined a numerical solution in the region where the 
gravitational curvature potential varies rapidly, with an approximate WKB solution for relatively large values of the 
radial coordinate. This trick allowed them to perform calculations for / < 20 and wM < 2.5. Some years ago one 
of us used the phase-integral method [TslJT^ to derive an approximate formula for the phase-shifts in the context of 
Schwarzschild scattering This formula was shown to be reliable and efficient even for high frequencies and/or 
large I values (in |^ results for luM — 10 and / < 200 were presented). This means that the differential cross sections 
determined from the phase-integral phase-shifts were reliable also for rather high frequencies. Even though the phase- 
integral formula could be generalised to scattering by a Kerr black hole and therefore used for the purposes of the 
present study, we have chosen a different approach here. Our phase-shift determination is based on the so-called Priifer 
method (well-known in quantum scattering theory l]!^,^ and, in general, in numerical treatments of Sturm-Liouville 
problems ||l9[| ) which, in a nutshell, involves transforming the original radial wavefunction to specific phase-functions 
and numerical integration of the resulting equations. In essence, this method is a close relative of the phase-amplitude 
method that was devised by one of us to study black-hole resonances |2^ . 

The remainder of the paper is organised as follows. In Sections IIA and IIB the problem of scattering by a Kerr 
black hole is rigorously formulated. In Section IIC the important notion of the deflection function is discussed. Section 
III is devoted to our numerical results. First, in Section IIIA our numerical method for calculating phase-shifts is 
presented. In Section IIIB familiar Schwarzschild results are reproduced as a code validation. Sections IIIC and HID 
contain entirely new information: Differential cross sections and deflection functions for on and off-axis scattering 
respectively. These are the main results of the paper. Furthermore, in Section HIE we present numerical results 
concerning forward glories. The role of superradiance for scattering of monochromatic waves is discussed in Section 
IIIF. Our conclusions are briefly summarised in Section IV. Three appendices are devoted to technical details, which 
are included for completeness. In Appendix A we discuss the notion of "plane waves" in the presence of a gravitational 
field. In Appendix B the partial-wave decomposition of a plane wave in the Kerr background is determined, and finally 
in Appendix C we briefly describe the method we have used to calculate the spin-0 spheroidal harmonics and their 
eigenvalues. Throughout the paper we adopt geometrised units (c = G = 1). 

II. SCATTERING FROM BLACK HOLES 

A. Formulation of the problem 

We consider a massless scalar field in the Kerr black-hole geometry. Then, first-order black-hole perturbation theory, 
basically the Teukolsky equation ||2^], applies. The scalar field satisfies the curved spacetime wave equation □$ = 0. 
Adopting standard Boyer-Lindquist coordinates we can always decompose the field as (since the spacetime is axially 
symmetric) 

^r,e,ip,t) = ^== J2 ^mir,e,t)e^"'^ (1) 
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In scattering problems it is customary to consider monocromatic waves with given frequency lo. Therefore we can 
further write 



+00 



(2) 



l — \m\ 

where cim is some expansion coefScient and S^^'^{9) are the usual spin-0 spheroidal harmonics. These are normalised 



as 



Finally, the function ui„i{r,uj) is a solution of the radial Teukolsky equation 

dPuirn , \ K'^ + {2amLo — a^u}"^ — Eirn)^ dG 
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where K = {r'^ + a^)uj — am and G ~ rA/{r'^ + a'^)'^. Furthermore, Eim denotes the angular eigenvalue, cf. Appendix C. 
As usual, A — r'^ — 2Mr + and the "tortoise" radial coordinate r, is defined as (with r±, the two solutions to 
A = 0, denoting the event horizon and the inner Cauchy horizon of the black hole) 



^^ln(- 



2Mr_ / r 
In 1 I + c 



(5) 



Usually, the arbitrary integration constant c is disregarded in this relation. However, in scattering problems it turns 
out to be useful to keep it, as we shall see later. 

We are interested in a causal solution to (Q) which describes waves that are purely "ingoing" at the black hole's 
horizon. This solution can be written 



as r 
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where k = u)~ma/2Mr^ = uj — muj^. In addition, we want to impose an "asymptotic scattering boundary condition". 
We want the total field at spatial infinity to be the sum of a plane wave plus an outgoing scattered wave. In other 
words, we should have 
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where we have omitted the trivial time-dependence. All information regarding scattering is contained in the (complex- 
valued) scattering amplitude f{9,Lp). Note that, unlike in axially symmetric scattering the scattering amplitude will 
depend on both angles: 9 and Lp. 

Up to this point, we have used the term "plane wave" quite loosely. In the presence of a long-range field such as 
the Kerr gravitational field (which falls off as ~ 1/r at infinity) we cannot write a plane wave in the familiar flat 
space form. This problem has been discussed in several papers, see Remarkably, it turns out that in a black 

hole background the long-range character of the field is accounted for by a logarithmic phase-modification of the flat 
space plane- wave expression. In practice, the substitution r — > is made in the various exponentials. In order to 
make this paper as self-contained as possible, we discuss this point in some detail in Appendix A. 



FIG. 1. A schematic illustration of the general scattering problem. A plane wave impinges on a rotating black hole making 
an angle 7 with the rotation axis. 
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The asymptotic expression for a plane wave travelling along a direction making an angle 7 with the black hole's 
spin axis, see Figure |l], is 

iojr* (sin 7 sin ^ sin c^+cos 7 cos 0) 

plane — yoj 

where without any loss ofgenerality we have assumed an amplitude of unity. We can decompose this plane wave in 
a way similar to (|^) and (||); 

+00 +00 

m=-oo l = \m\ 

where are asymptotic solutions of (Q). For r oo we have (see Appendix B), 

cXir,^) « '2nSr-h) +^'"+l(„l)'+'"e— } (10) 

Similarly, the full field at infinity can be approximated as : 



+00 +00 
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By imposing the scattering condition (^, we can fix by demanding that the ingoing wave piece of $ — $piane 
vanishes. After some straightforward manipulations we get for the scattering amplitude 



„ +00 +CXD 
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By defining the "scattering matrix element" Sim = (^l)'^^^°m*/^/m '^^ can equivalently write 

Sim = e^^*'" (13) 

where we have introduced the phase-shift 5im- Thus we see that the phase-shifts 5im contain all relevant information 
regarding the scattered wave. It is worth emphasising that for non-axisymmetric scattering the phase-shifts will 
depend on both / and m. Also, the 5im are in general complex valued in order to account for absorption by the black 
hole. For later convenience, we also point out that the full field at infinity can be written 

$ ~ sin(ijjr, + 5im — —) as — + cxd (14) 
In the case of on-axis incidence (7 = 0) the scattering amplitude simplifies considerably, and we get 

(15) 



/w = -E^"o"mT(o) 
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Here we see that the outcome is no longer dependent on m, which is natural given the axial symmetry of the problem. 
Furthermore, it is easy to see that we recover the familiar Schwarzschild expression ^ by setting a = 0. 

The difi'erential cross section (often simply called the cross section in this paper) is the most important "observable" 
in a scattering problem. It provides a measure of the extent to which the scattering target is "visible" from a certain 
viewing angle. As demonstrated in standard textbooks [0, the differential cross section follows immediately from 
the scattering amplitude 

^ = i/(e,¥')r (16) 

This cross section corresponds to "elastic" scattering only, that is, it describes the angular distribution of the waves 
escaping to infinity. We can similarly define an "absorption cross section" but we shall not be concerned with this 
issue here. Nevertheless, as we have already pointed out, black hole absorption has an effect on the phase-shifts that 
are used to compute the cross section (p^. 

The strategy then for a cross section calculation (for given black hole parameters and wave frequency) involves three 

steps: i) calculation of the phase-shifts Sim, (or, equivalently, of the asymptotic amplitudes ii) calculation of 

the spheroidal harmonics 8^(9), and finally iii) evaluation of the sums in ( p^ ) and/or ( p^ ) including a sufficiently 
large number of terms. 
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B. Approximating the scattering amplitude 



In practice, the partial-wave sum calculation is problematic as it converges slowly. In fact, the sum is divergent for 
some angles. This is just an artifact due to the long-range nature of the gravitational field. No matter how far from 
the black hole a partial wave may travel, it will always "feel" the presence of the gravitational potential (that falls 
off as l/r) . A similar behaviour is known to exist in Coulomb scattering. The divergence always occurs at the angle 
that specifies the incident wave's propagation direction. That this will be the case is easily seen from the identity 
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(7)e" 
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5{cos6 — cos j)S{(p — tt/2) 



(17) 



which follows directly from the fact that the functions 5) 
identity for on axis incidence is 



(6')e™'^ form an orthonormal set. The corresponding 
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mT(O) = ^<5(cos0-l) 



(18) 



From ( [l7| ) we can deduce a peculiar feature: Although the scattering problem is physically insensitive to the actual 
if of the incidence direction, the specific value ip — 7r/2 is imposed by the above relations. Of course, this has no 
physical relevance since the problem at hand is axially symmetric and we can, without any loss of generality, assume 
an incoming wave travelling along the direction {0,ip) — (7,7r/2). 

The fact that the Kerr gravitational field behaves asymptotically as a Newtonian one considerably simplifies the 
scattering amplitude calculation. We would expect that large I partial waves (strictly speaking when I/loM ^ 1) 
to essentially feel only the far-zone Newtonian field. In terms of the phase-shifts, we expect them to approach their 
Newtonian counterparts Sim — ^/^ asymptotically. In order to secure this matching we add to our phase-shifts an 
"integration constant" —2u)Mhi{'^LoM) + uiM. In this way, we also get — > Tc, where rc = r + 2Mln(2wr) is 
the respective tortoise coordinate of the Coulomb/Newtonian problem. Such a manipulation is admissible given the 
arbitrariness in the choice of the constant c in (^ . 

In calculating the partial-wave sum for the scattering amplitude, it is convenient to split it into two terms: 



f{0,cp) = fD{0,v)+fN{O,v) 



(19) 



Here, foid,^) represents the part of the scattering amplitude that carries the information of the main diffraction 
effects, while fN{0,(p) denotes the Newtonian (Coulomb) amplitude. Explicitly we have 



27r 



^l^™(e)e""^r/™(7)(-^)^ 
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(20) 



where we have deliberately "forgotten" the spherical symmetry of the Newtonian potential (which would had allowed 
us to write Jn as a function of 6 only, and thus in terms of a sum over /). However, the Newtonian phase-shifts Sj^ 
are still given by the well-known expression [||, 



After simple manipulations we get 
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where cos^ — cos 6' cos 7 + sin sin 7 sin The sum in (g2|) is known in closed form ||1[|; 



r(l - 2iLuM) 
r(l + 2iujM) 
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(23) 



From this we can see that /Ar(6', ip) diverges in the ^ — direction. 

Let us now focus on the "diffraction" amplitude f]j{0,ip). It has the form 
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The corresponding on-axis expression is, 

+00 

/^(^) = ^ E {47r5ro" W5ro"(0)(e2''' - 1) - [21 + mico,e)ie^< - 1)} (25) 



One would expect the sums in (24) and (ES) to converge. This foUows from the fact that for l/ujM ^ cx) we have 
5im — > if and Sf^{9)e™'"^ Yim{0,ip) [||~^We introduce a neghgible error by truncating the sums at a large value 
^max (say). In practice, /max need not be very large. We find that a value ~ 30 — 50 for loM < 2 typically suffices. 
Since each partial wave can be labelled by an impact parameter b{l) (see Section IID), the criterion for /max to be a 
"good" choice, is that 6(/niax) S> 60, where 6c is the (largest) critical impact parameter associated with an unstable 
photon orbit in the Kerr geometry. 

The truncation of the partial-wave sums will introduce interference oscillations in the final cross sections (roughly 
with a wavelength 27r//max (^])- These unphysical oscillations can be eliminated by following the approach of Handler 
and Matzner 0. For a chosen l^ax we add a constant (3 to all the phase-shifts in (p3) and (|2^). This constant is 



chosen such that (5;,„^^,m + 13 = Sf^ . This means that the resulting cross section is effectively smoothed. 



C. Deflection functions 



It is well- know that the so-called deflection function is of prime importance in scattering problems. It arises in the 
semiclassical description of scattering, as discussed in the pioneering work of Ford and Wheeler . Although these 
authors considered scattering in the context of quantum theory, their formalism is readily extended to the black-hole 
case. In the semiclassical paradigm, the phase-shifts are approximated by a one-turning point WKB formula (typically 
useful for I much larger than unity). 

In a problem which has only one classical turning point, the deflection function is deflned as 



6(0 



' dl 



(26) 



where I is assumed to take on continous real values. As a convention, the deflection function is negative for attractive 
potentials. The right-hand side of this equation resembles the expression for the deflection angle of classical motion 
in the given potential, provided that we define the following effective impact parameter b for the wave motion pi] 



/ + 1/2 



(27) 



The black-hole effective potential has two turning points, but for 1/ljM 00 the scattering is mainly due to 
the outer turning point and one can derive a one turning point WKB approximation for the phase shifts. For a 
Schwarzschild black hole this expression is |q| 



+ 00 



Qs 



_/ _2MY 



r / 



dr ~ IjjU + {21 + 1)^ 
4 



(28) 



where 



Ql 



2M 
r 



2M 



£\ 1(1 + I) 



r J 



(29) 



Here denote the value of the tortoise coordinate corresponding to the (outer) turning point t. We now define the 
deflection function as 



e(/) = 2Re 



dS, 



dl 



(30) 



where only the real part of the phase shift is considered, as the whole discussion is relevant for elastic scattering only. 
Using (p8|), we find that the real scattering angle is 
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9(0 



1 + 1/2 



2M\ f 1 + 1/2 



-0{NP/r^ 



-1/2 



(31) 



This WKB result should be compared to the deflection angle for a null geodesic in the Schwarzschild geometry, which 
is given by 



dr 

ee(6)=7r-26/ ^ 



1-1- 



2M\ h 



1-1/2 



(32) 



where b = L^/E is the orbit's impact parameter (L^ and E denote, respectively, the orbital angular momentum 
component along the black hole's spin axis and the orbital energy) and t is the (classical) turning point. In writing 
down these expressions we have chosen the signs in such a way that the deflection angle is negative for attractive 
potentials. Clearly, it is possible to "match" the deflection function (^) with the classical deflection angle, albeit 
only at large distances. Since the effective impact parameter will be given by (|2^), it is clear that in (^) the integral 
will be over large r only. 

Owing to its clear geometrical meaning the deflection function is an exceptionally useful tool in scattering theory. It 
can be used to define diffraction phenomena like glories, rainbows etc. |25|] . For example, in axisymmetric scattering 
backward glories are present if the deflection function takes on any of the values O = — titt, where n a positive odd 
integer. 

It seems natural to try and define deflection functions for Kerr scattering as well. In general, we anticipate the 
need for two deflection functions 0(^, m) and ^{l, m) (with only the first being relevant for the special case of on-axis 
scattering). The WKB phase-shift formula becomes in the Kerr case: 



rWKB 



+ 00 



Qk- 



dr - Lot^ + (21 + 1)- 
4 



where 



A2 



(33) 



(34) 



where A — Eim + a^uj'^ — 2amuj. 

The next step is to derive the deflection angles for null geodesies approaching a Kerr black hole from infinity. Such 
orbits are studied in detail in j2^. From the results in j2^ it is clear that it is not easy to write down a general 
expression for the defiection angle in the Kerr case. But we can obtain useful results in two particular cases. 

We begin by considering a null ray with = (which would correspond to an axially incident partial wave) . For 
such an orbit we find that the deflection angle 8c obeys the following relation 



de 



1 + ^ cos^ 



-1/2 



. dr 
-2^ Z7 



1-^ 



1 



2M 
r 



-1/2 



(35) 



where ry = C^^'^/E, with C denoting the orbit's Carter constant. This expression is valid provided the ray's 9- 
coordinate varies monotonically during scattering. This should be true in the cases we are interested in, at least for 
large impact parameters such that 77 ^ M. The ray will also be deflected in the (^-direction but this deflection carries 
no information regarding plane-wave scattering due to the axisymmetry of the problem. 

We next consider a null ray travelling in the black hole's equatorial plane. This situation will be particularly relevant 
for a plane wave incident along 7 = 7r/2. The net azimuthal deflection <I'c(&) for impact parameter b — L^/E is 



$c(6) =^-26^ 



dr 
r 



2aMr 

A ^ 



A5 



. , 2a'^M 2M\ 4aMb 
+a^ + feM 1 - — j 



1-1/2 



(36) 



Working to the same accuracy in terms of M/r as in the Schwarzschild case, we can match ( ^5|) and ( ^) to d5im/dl. 
This matching becomes possible if we use the following approximate expression for the eigenvalue M 



Ei^^l{l + l)- -o'lo^ + 0{- 



l 



(37) 



As in the Schwarzschild case we assume that the effective impact parameter is given by (^^. Although there is no 
occurrence of the multipole m in the above expressions, one can argue (from the symmetry of the various spheroidal 
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harmonics, which is similar to that of the spherical harmonic of the same {I, m)) that the classical angles ( ^5| ) and ( |3^ ) 
are related to partial waves with to = and incidence 7 = and partial waves with m = ±1 and incidence 7 — 7r/2, 
respectively. Hence, we define the latitudinal deflection function 



6(0 = 2Re 



and the azimuthal ("equatorial") deflection function 



dl 



(to = 0) 



$(0 = 2Re 



dSi, 



dl 



-(m = ±1) 



(38) 



(39) 



III. NUMERICAL RESULTS 



A. Phase-shifts calculation via the Priifer transformation 



In order to determine the required scattering phase-shifts we have used a slightly modified version of the simple 
Priifer transformation, well-known from the numerical analysis of Sturm-Liouville problems p9[ | . The method is best 
illustrated by a standard second order ordinary differential equation: 



(40) 



where we can think of x as being a radial coordinate, spanning the entire real axis, and U an effective potential (in 
our problem corresponding to a single potential barrier) with asymptotic behaviour 



Uix) 



as X - 
as X 



'OO 

I- 00 



(41) 



with to and k real constants. (The black hole problem we are interested in does, of course, have exactly this nature.) 
The solution of (Eoh will take the form of oscillating exponentials for x ±00. Let us assume that we are looking 
for a solution to (HO) with purely "ingoing" behaviour at the "left" boundary {x — > —00) and mixed ingoing/outgoing 
behaviour as x +co: 



e ''^^ as X 



-00 



B siii[Lux + C] as X — > +00 . 
where ^ and B are complex constants. We can then write the exact solution of (^0|) in the form 

The function P{x) is the logarithmic derivative of ip{x) (a prime denotes derivative with respect to x) 

V'' 

— = p 

which obeys the boundary condition P{x) —ik for x — > —00. 

Similarly, we can express the function tp and its derivative via a Priifer transformation; 



ip{x) ~ Bsm[ujx + P{x)] 
ip'{x) = Buj cos[llJX + P{x)] 

with P{x) a Priifer phase function which has C as its limiting value for x 
the equations 



(42) 



(43) 



(44) 



(45) 
(46) 



^ + P^ + U{x) 
ax 







-00. Direct substitution in ( pO| ) yields 

(47) 
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dP 

dx 



U{x) 



ain^ [P + ujx) = 



(48) 



The idea is to numerically integrate (^) and (^8|) instead of the original equation (|4^). The motivation for this is 
that, while the original solution may be rapidly oscillating, the phase-functions P and P are expected to be slowly 
varying functions of x. 

We expect this integration scheme to be considerably more stable, especially for high frequencies, than any direct 
approach to (40). Moreover, eqs. (|4^) and (|4^) are well behaved also at the classical turning points and are well 
suited for barrier penetration problems. However, if we want to ensure that the phase-functions are smooth and 
non-oscillatory we must account for the so-called Stokes phenomenon — the switching on of small exponentials in 
the solution to an equation of form (^). To do this we simply shift from studying P{x) (which is calculated from 
X — — oo up to the relevant matching point Xm) to P{x) (which is calculated outwards to a; = +cxd). In practice, 
the calculation is stable and reliable if the switch is done in the vicinity of the maximum of the black-hole potential 
barrier (the essential key is to not use one single representation of the solution through the entire potential barrier). 
The two phase functions are easily connected by 



P{x) 



-UJX - 



2i 



■In 



iP 



iP 



(49) 



Finally, the desired phase-shift can be easily extracted as 6im = C + Itt/2. 

A major advantage of the adopted method is that it permits direct calculation of the partial derivatives dSim/dl 
and dSim/dm, which are required for the evaluation of the deflection functions from Section IID. The equations for the 

m-derivatives of P and P are simply found by differentiation of (47) and (^8|). For the case of scattering by a Kerr 
black hole, calculation of 5im and its derivatives with respect to I, m requires knowledge of the angular eigenvalue Eim 
(see Appendix C) and its_J, m-derivatives. We have used an approximate formula which is a polynomial expansion 
in auj (formula 21.7.5 of |^^). This expression (and its derivatives) is well behaved for all integer values of I and to. 
However, it is divergent for the half-integer values I = 1/2, 3/2, 5/2. Hence, the numerical calculation of the deflection 
function will fail at these points, and will be generally ill-behaved in their neighbourhood. For the full cross section 
calculation, we additionally need to calculate the spin-0 spheroidal harmonics Sf^{9). This calculation is discussed 
in detail in Appendix C. 



B. Schwarzschild results 



In this section we reproduce phase-shifts and cross sections for Schwarzschild scattering. The purpose of this exercise 
is to validate, and demonstrate the reliability of, our numerical methods. We compare our numerical integration results 
to ones obtained using the phase-integral method [||. 

As a first crucial test we compare, in Fig. ||, the first 100 phase-shifts for u)M — 1. Because of the multi-valued 
nature of the phase-shifts we always plot the quantity Si = e^*''' . As is evident from Fig. ||, the agreement between 
our numerical phase-shifts and the phase-integral ones is excellent. This is equally true for a all frequencies examined 
(up to LoM = 10). It should be noted that Si is essentially zero {5i has a positive imaginary part) for those values 
of / for which absorption by the black hole is important. That this is the case for the lowest multipoles is clear from 
Fig. |. As I increases 5i becomes almost real and as a consequence Si is almost purely oscillating. 
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FIG. 2. Comparison of numerical phase-shifts for a Schwarzschild black hole (cross), against phase-integral data (plus). We 
show the real (upper frame) and imaginary (lower frame) parts of the scattering matrix element Si = e^'*' as functions of 1. 
The agreement between the two methods is clearly excellent. 

In Fig. ^ we present an ojM — 10 cross section generated from our numerical phase-shifts. For this particular 
calculation we have used /max — 200. The resulting cross section matches the one constructed using approximate 
phase-integral phase shifts perfectly. This demonstrates the efficiency of our approach in the high frequency regime, 
and it is clear that our study of Kerr scattering will not be limited by the lack of reliable phase shifts. However, 
the Kerr study is nevertheless limited in the sense that uj cannot be taken to be arbitrarily large. This restriction is 
imposed by the calculation of the spheroidal harmonics (Appendix C). However, it is important to emphasize that 
the most interesting frequency range, as far as diffraction phenomena is concerned, is ujM 2 Thus, we expect 
that our investigation should be able to reliable unveil all relevant rotational effects in the scattering problem. 
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FIG. 3. Differential cross section for scattering of a wave with relatively high frequency, coM — 10, from a Schwarzschild 
black hole, based on the the first 200 partial wave phase-shifts. The backward glory oscillations are prominent. 

We also find that the numerical values for d6i/dl are in good agreement with the respective phase- integral results. 
As a final remark we should emphasize that the numerical approach adopted in this work is very efficient from a 
computational point of view. 



C. On-axis Kerr scattering 



Having confirmed the reliability of our numerical results we now turn to the study of scattering of axially incident 
scalar waves by a Kerr black hole, cf. Fig. |4[ In principle, we would expect the corresponding cross sections to be 
qualitatively similar to the Schwarzschild ones. The main reason for this is the inability of axially impinging partial 
waves to distinguish between prograde and retrograde orbits. However, examination of the orbital equations p^ ] 
reveals that the critical impact parameter (associated with the unstable photon orbit) decrease sligthly from the value 
SVSAf as the black hole spins up. For example, for a=0.99M we have be = 4.74M. 
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FIG. 4. A schematic drawing illustrating the case of on-axis scattering from a rotating black hole. Because of the axial 
symmetry of the problem, the scattering results are qualitatively similar to those for a Schwarzschild black hole. Two rays, 
which emerge having been scattered by the same angle are indicated. 

Indeed, Fig. ^ confirms our expectations. The data in the figure corresponds to a black hole with spin a = 0.99M 
and a wave frequency of loM = 2. For purposes of comparison, we also show the corresponding Schwarzschild cross 
section. The two cross sections are very similar. In particular, they are both dominated by the backward glory. It 
is well-known |^J^ that for Schwarzschild scattering the glory effect can be described in terms of Bessel functions. It 
has been shown that for 9 n the glory cross section can be approximated by 



^(6*) I glory oc Jo[w6cSin6'] 



(50) 



A similar result holds for the backward glory in the case of on-axis Kerr scattering. Since be gets smaller, one would 
expect the zeros of the Bessel function (the diffraction minima) to move further away from 9 — tt as a increases. This 
effect is indicated by the data in Fig. 
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FIG. 5. Differential cross section for on-axis scattering from a Kerr black hole (solid line). The black hole's spin is a = 0.99Af 
and uiM — 2. The graph is based on data for Imax ~ 30. The dashed line represents the corresponding Schwarzschild cross 
section. 



The above conclusions are further supported by the results for the deflection function (^8|), as shown in Fig. ||. 
Because of the inaccuracies inherent in our method of calculating the deflection function for the lowest Z-multipoles, see 
the discussion in Section IIIA, we do not show results for this regime. This is, however, irrelevant as the corresponding 
partial waves are expected to be more or less completely absorbed by the black hole. 
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FIG. 6. Deflection function (in units of tt) for on-axis scattering from a Kerr black hole. The black hole spin is a = 0.99M 
and uM = 1. The dashed curve is the corresponding Schwarzschild deflection function. This figure confirms the behaviour 
expected from the geometric optics considerations, namely, the slight decrease of the critical impact parameter with increasing 
a. 



D. Off-axis Kerr scattering 

The conclusions of our study of on-axis scattering are perhaps not very exciting. Once the Schwarzschild case is 
understood, the on-axis results for Kerr come as no surprise. This is, however, not the case for off-axis scattering, cf. 
Fig. 1^, where several new features appear. 




(p=0 

FIG. 7. A schematic drawing illustrating off-axis scattering from a rotating black hole. We show (as thick dashed lines) two 
"rays", one of which corresponds to motion in the black hole's equatorial plane. 

Our study of the off- axis case provides the first results for non-axisymmetric wave scattering in black hole physics. 
Since this problem has not been discussed in great detail previously, it is worthwhile asking whether we can make any 
predictions before turning to the numerical calculations. Two effects ought to be relevant: First of all, the partial 
waves now have orbital angular momentum which couples to the black hole's spin. As a result the partial waves can 
be divided into prograde (m > 0) and retrograde (m < 0) ones. We expect prograde waves to be able to approach 
closer to the horizon than retrograde ones. In the geometric optics limit, prograde and retrograde rays tend to have 
increasingly different critical impact parameters as a — ^ M. As a second feature, we expect to find that large I partial 
waves will effectively feel only the spherically symmetric (Newtonian) gravitational potential. In other words, partial 
waves with the same (large) I and different values of m will approximately acquire the same phase-shift. 

Our numerical results essentially confirm these expectations, as is clear from the phase-shifts (calculated for a = 
0.9M and ujM = 1) shown in Fig. ^. As above, we have graphed the single-valued quantity Sim = e^*''''" as a function 
of I. For each value of I we have included all the phase-shifts for —l<m< +1. The solid (dashed) line corresponds 
to m = +^ (to = —I) and the intermediate values of to lead to results in between these two extremes. For I 3> 1, 
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partial waves with different values of m have almost the same phase-shift. This is easy to deduce from the fact that 
the two curves approach each other as / increases. On the other hand, for the first ten or so partial waves we get 
very different results for the various values of m. In particular, we see that phase-shifts with m > become almost 
real (that is, \Sim\ becomes non-zero) for a smaller Z-value as compared to the m < ones. As anticipated, this is 
due to the different critical impact parameters associated with prograde/retrograde motion, and the fact that a larger 
number of prograde partial waves are absorbed by the black hole. 




We now turn to the cross section results for the off-axis case. We have considered a plane wave incident along the 
direction ^ — (p = 7r/2. Even though our formalism allows incidence from any direction we have focussed on this case, 
which is illustrated in Fig. ^. The motivation for this is that there will then be partial waves (specifically the ones with 
m = ±/) that are mainly travelling in the black hole's equatorial plane. These partial waves are important because 
one would expect them to experience the strongest rotational effects. Besides, we can obtain an understanding of 
these waves by studying equatorial null geodesies in the geometric optics limit. Equatorial null rays are much easier 
to describe than nonequatorial ones. This proves valuable in attempts to "decipher" the off-axis cross sections, and 
the obtained conclusions provide an understanding also of the general case. 

In Fig. ^ we present a series of cross sections as functions of ip for the specific values 6 — 7r/8, 7r/4, 37r/8, 7r/2. 
These results correspond to viewing the scattered wave on the circumference of cones (like that shown in Fig. ^ with 
increasing opening angles. Two different frequencies ujM = 1 and uiM = 2 have been considered for a black hole with 
spin a — 0.9Af . A first general remark concerns the asymmetry of the cross sections with respect to the incidence 
direction (note, however, that as a consequence of our particular choice of incidence direction there is still a reflection 
symmetry with respect to the equator) . We can also easily distinguish the Coulomb forward divergence in the direction 
= Lp = i: /2. Another obvious feature in Fig. ^ is the markedly different appearance of the cross section for different 
values of 9. As we move away from the equatorial plane the cross sections becomes increasingly featureless. This 
behaviour is artificial in the sense that as 9 decreases, we effectively observe along a smaller circumference. Kl 6 — 
this circumference degenerates into a point, cf. Fig. 0. 
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FIG. 9. Off-axis cross sections {6 — vr/S, 7r/4, Stt/S, 7r/2) for a black hole with spin a = 0.9M and scattered waves with 
frequency uiM = 1 (upper panel) and uM — 2 (lower panel). The incident wave is travelling in the {0,if) = (7r/2,7r/2) 
direction. 

In order to understand the features seen in Fig. ^further, we focus on the 6 = tt/2 cross section. In Fig. ^we show 
these "equatorial" cross sections for a sequence of spin rates a/M — 0.2,0.5,0.7,0.9. As before, we have considered 
two different wave frequencies, coM = 1 and ujM = 2. From the resuhs shown in Fig. |l^ it is clear that that the glory 
maximum is typically not observed in the backward {(p = —tt/2) direction. In fact, it is clear that the maximum of 
the glory oscillations move away from the backward direction as the spin of the black hole is increased. A similar 
shift is seen in all interference oscillations. This behaviour is easy to explain in terms of the anticipated rotational 
frame-dragging. In order to illustrate this argument, we consider the geometric optics limit where partial waves are 
represented by null rays. Recall that in axisymmetric scattering the backward glory is associated with the divergence 
of the classical cross section at 9 = tt in such a way that 

\dn sine \ db ) ^ ' 



The divergence is a result of the intersection of an infinite number of rays. For simplicity, let us consider rays travelling 
in a specified plane. Scattering near the backward direction by a Schwarzschild black hole is illustrated in Fig. O 



15 



(left panel). Two rays with different impact parameters emerge at any given angle. These two waves make the main 
contribution to cross section at that particular angle. It is clear that for scattering a,t 6 — tt the two rays in Fig. |ll| 
will follow symmetric trajectories. This means that when observed at infinity, after being scattered, the two waves 
will have equal phases (provided their initial phases were equal). In effect, these two rays will then constructively 
interfere in the exact backward direction. As we move away from the backward direction, we should observe a series of 
interference maxima and minima — the two rays will now have an overall phase difference since they follow different 
orbits (see Fig. 11). A very crude estimate of the location of the successive maxima would be 9^ ^ mr/SujM where 
n = 0, 1, 2, in reasonable agreement with the exact results. 





-0.5 00 

(p/7t 

FIG. 10. Off-axis cross sections for 9 — n/2 and various black-hole spins {a,/M— 0.2, 0.5, 0.7, 0.9) and scattered wave 
frequencies uM = 1 (upper panel) and ujM = 2 (lower panel). 

Similar arguments apply in the case of a Kerr black hole. We shall consider only equatorial rays, cf. Fig |ll| (right 
panel). As a result of the discrimination between prograde and retrograde orbits, the two rays contributing to the 
cross section in the exact backward direction will no longer follow symmetric paths. In fact, the ray symmetric to the 
prograde ray shown in Fig[l^ will follow a plunging orbit. Therefore, we should not expect the interference maximum 
to be located in the exact backward direction. An estimate (based on a crude calculation of the phase difference 
between the two null rays) of the location of the main backward glory maximum yields 
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^max - n (!Ph±^:i!hA (52) 

where rph+ and fph- denote, respectively, the location of the prograde and retrograde unstable photon orbits (in 
Boyer-Lindquist coordinates). This angle is measured from the backward direction in the direction of the black hole's 
rotation. This simple prediction agrees reasonably well with the results inferred from our numerical cross sections. In 
a similar way, all other maxima and minima will be frame-dragged in the black hole's rotational direction. 
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FIG. 11. Equatorial null geodesies (viewed from "above") around a Schwarzschild (left panel) and Kerr black hole (right 

panel). The figures are scaled in units of M. The rays are assumed to arrive from infinity (they enter from the right side of 

each figure in the direction indicated by the arrows) in parallel directions and exit at the same angle after being scattered. 

The dashed circles represent the unstable photon circular orbits. The Kerr black hole, in the right panel, is taken to rotate 

counter-clockwise with a = 0.99M. 

To complete this discussion we consider the deflection function <^(l,m) for "equatorial" partial waves (m = ±^), an 
a — 0.9M black hole and ivAI = 1. The corresponding data is shown in Fig. n2[ There are two distinct logarithmic 
divergences which are associated with the existence of separate unstable circular photon orbits for prograde and 
retrograde motion. Note that for m > the deflection function diverges steeper than it does for m < 0. The origin 
of this effect is the fact that prograde partial waves with b ^ be perform a greater number of revolutions (before 
escaping to infinity) than retrograde ones. Finally, for |m| ^ 1 we recover, as expected, the Einstein deflection angle 
$ « —AM/b (not explicitly shown in the figure). 
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FIG. 12. The defiection function $ for off-axis scattering is shown as a function of m for "equatorial" partial waves m — ±L 
The black hole spin is a = 0.9M and the wave frequency is luA-I — 1. The small |m| region is not included because of the 
inaccuracies discussed in Section III A. 



E. Digression: forward glories 

The results presented in the preceeding sections clearly show that, in general, black hole cross sections are dominated 
by a "Coulomb divergence" in the forward direction and (frame-dragged) glory oscillations near the backward direction. 
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However, according to the predictions of geometrical optics ||l[, one would expect to find glory oscillations also in the 
forward direction (see comments in |^). For the case of Schwarzschild scattering, this effect would be associated with 
partial waves scattered at angles 9 = 0, — 27r, —An, .... Inspection of the relevant deflection function (Fig. ^ indicates 
that a partial wave which has = will also be strongly absorbed, since it has an impact parameter h < h^. Hence, 
we would expect its contribution to the forward glory to be severely supressed. It thus follows that, as far as the 
possible forward glory is concerned, the most important partial waves are those with = — 27r. These partial waves 
"whirl" around the black hole as they have h k, be. 

Ford and Wheeler's semiclassical approach shows that the forward glory is well approximated (for 6* w 0) 
by (50), although with a slightly different proportionality factor. However, we should obviously not expect to see 
a pronounced forward glory in the cross section, as it will drown in the forward Coulomb divergence. Still, as an 
experiment aimed at supporting our intuition, we can try to "dig out" the forward glory pattern. This has to be done 
in a somewhat artificial manner, but since the forward glory is due to scattering and interference of partial waves 
with h ~ be we can isolate their contribution by truncating the sum in at some Zmax ~ w&c and at the same time 
neglecting the Newtonian part /at entirely. It is, of course, important to realize that this "truncated" cross section is 
not a physical (observable) quantity. 

In Fig. |l^ we show the result of this "truncated cross section" calculation for the case of a Schwarzschild black 
hole and loM — 2. We compare results for two levels of truncation, /,„ax = 10 and 15. In the first case, a clear 
Bessel- function like behaviour arises (it is straightforward to fit a Jq (cjfec sin 6*) function to the solid curve in Fig. p^ . 
This confirms our expectation that there is, indeed, a forward glory present in the data. As more partial waves are 
included the cross section begins to deviate from the glory behaviour, and if we increase ^,nax further the forward 
glory is swamped by terms that contribute to the Coulomb divergence. 
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FIG. 13. Illustration of a forward glory in Schwarzschild scattering. The "diffraction" piece |/D(e)|2 of the cross section is 
shown in the vicinity of the forward direction, for wave frequency ujM — 2 and for /max = 10 (solid curve) and /max = 15 
(dashed curve). 



F. The role of superradiance in the scattering of monochromatic waves. 



Superradiance is an interesting effect known to be relevant for rotating black holes. It is easily understood from the 
asymptotic behaviour (^ of the causal solution to the scalar- field Teukolsky equation (^) . If we use this solution and 
its complex conjugate, and the fact that two linearly independent solutions to (j^ must lead to a constant Wronskian, 
it is not difficult to show that 

(l-mco+M|7J„.|2 = 1- |5,„J2 . (53) 

where we have defined 

(54) 
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From the above result it is evident that the scattered waves are ampHfied > 1) if w < mu}+. This amphfication 

is known as superradiance. 

In principle, one would expect superradiance to play an important role in the scattering problem for rapidly 
spinning black holes. For example, one could imagine that some partial waves which would otherwise be absorbed, 
could escape back to infinity. These waves might then possibly make a noticeable contribution to the diffraction cross 
section, provided that there were a sufficient number of them (as compared to the total number ^max of partial waves 
contributing to the diffraction scattering amplitude). 

In order to investigate this possibility, we have performed a number of off-axis cross section calculations for a 
variety of wave frequencies ojM = 0.5 — 10 and for a « M, i.e. black holes spinning near the extreme Kerr limit. We 
have found no qualitive difference whatsoever between those cross-sections and the ones for a somewhat smaller spin 
value, a = 0.9M (say). In essence, we were unable to find any effects in the cross section that could be attributed to 
superradiance. Consequently, we are led to suspect that our intuition regarding the importance of superradiance for 
the scattering problem may be wrong. 

This suspicion is confirmed by the following simple argument. In order for a partial wave to be superradiant we 
should have Q < lo < muj+. Considering an extreme Kerr black hole (which provides the best case for superradiant 
scattering) and the fact that m < I, we have the condition 

< 2ujM < I (55) 

As already mentioned, the partial waves for which superradiance will be important are the ones with impact parameters 
b < be, i.e. those that would be absorbed under different circumstances. This then requires that 

1 <ujbe- 1/2 (56) 



Combining (55) and ( |56D we arrive at the inequality 

< 2M < 6c - l/2u; (57) 



The critical impact parameter (for prograde motion) for an a = M black hole is be ~ 2M. Hence the condition (|57_ 
will not be satisfied, and it is unlikely that we would get a significant number of (if any) superradiant partial waves 
that could affect the cross section. 

This conclusion may seem surprising given results present in the literature 0,^. In particular. Handler and Matzner 
have briefly discussed the effect of superradiance on axially incident gravitational waves. They argue that (see figure 14 
in 1^ ) "superradiance has the effect of imposing a large background over the pattern, filling in the interference minima" . 
Given our current level of understanding (or lack thereof) we cannot at this point say whether superradiance can be 
the explanation for the effects observed by Handler and Matzner. After all, one should remember that superradiant 
scattering strongly depends on the spin of the field that is being scattered. It is well known ||2^ that gravitational 
perturbations can be amplified up to 138% compared to a tiny 0.04% amplification for scalar fields (which is the case 
considered in this paper). This means that superradiance may significantly affect also partial waves with b > be in 
the gravitational wave case, which could lead to our simple argument not being valid. This issue should be addressed 
by a detailed study of the scattering of gravitational waves from rotating black holes. 



IV. CONCLUDING DISCUSSION 



We have presented an investigation of scattering of massless scalar waves by a Kerr black hole. Our numerical 
work is based on phase-shifts obtained via integration of the relevant radial wavefunction with the help of the Priifer 
phase-function method. This method has been shown to be computationally efficient and to provide accurate results. 
Using the obtained phase-shifts we have constructed differential cross sections for several different cases. First we 
have discussed the case of waves incident along the black hole's rotation axis, for which we showed that the resulting 
cross sections are similar to ones obtained in the (non-rotating) Schwarzschild case. We then turned to the case of 
off-axis incidence, where the situation was shown to change considerably. In that case the cross sections are generically 
asymmetric with respect to the incidence direction. The overall diffraction pattern is "frame dragged" , and as a result 
the backward glory maximum is shifted along with of the black hole's rotation. Moreover, we have concluded that (at 
least for scalar waves) the so-called superradiance effect is unimportant for monochromatic scalar wave scattering. 

To summarize, our study provides a complete understanding of the purely rotational effects involved in black-hole 
scattering. Given this we are now well equipped to proceed to problems of greater astrophysical interest, particularly 
ones concerning gravitational waves. In these problems one would expect further features to arise as the spin and 
polarisation of the impinging waves interact with the spin of the black hole. For incidence along the hole's spin axis. 
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one can have circularly polarised waves which are either co- or counter-rotating. The two cases can lead to quite 
different results. Although the general features of the corresponding cross sections are similar, they show different 
structure in the backward direction [|l|. This is possibly due to interference between the two polarisation states 
of gravitational waves, an effect that has not yet been explored in detail. Some initial work on gravitational-wave 
scattering has been done, see but we believe that the results of the present paper sheds new light on previous 
results, and could help interpret the rather complex cross sections that have been calculated in the gravitational-wave 
case. 

In this context, it should be stressed that the choice of studying scalar waves was made solely on grounds of clarity 
and simplicity. Our approach can readily be extended to other cases. Moreover, it is relevant to point out that a full 
off-axis gravitational wave scattering cross section calculation is still missing. We would expect such cross sections 
to be rather complicated, combining the frame-dragging effects discussed in this paper with various spin-induced 
features. We hope to be able to study this interesting problem in the near future. 
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APPENDIX A: PLANE WAVES IN THE KERR GEOMETRY 



The long-range character of the gravitational field modifies the form of "plane waves" . This non-trivial issue has 
been discussed in the context of black hole scattering by Matzner and Chrzanowski et al. |2^. For completeness, 
we provide a brief discussion here. 

In a field-free region a monochromatic plane wave is, of course, given by the familiar expression 



plane 



Aujr cos 6 — iujt 



(Al) 



when a spherical coordinate frame is employed. The plane wave is taken to travel along the z-axis. The field in (Al) 
solves the the wave equation Deplane = 0. Moreover, we can assume a decomposition of the form 



'i'plane 

where the radial wavefunction satisfies 



2 Ki + i) 



,(0) 



(A2) 



(A3) 



Let us now consider a "plane wave" in the Schwarzschild geometry. First of all, we expect such a field to be only 
an asymptotic solution (as r — > oo) of the full wave equation □$ — (where □ represents the covariant d'Alembert 
operator). The plane wave field can then be represented at infinity as 



plane 



■iujt 



(A4) 



and the radial wavefunction will be a solution of 



O 



In 



,(0) 



(A5) 



This equation is similar to the corresponding fiat space equation. The only difference is the appearance of the tortoise 
coordinate r* instead of r. Hence, we are inspired to write the plane wave field as 



plane 



^iuir^ cos 9 — iujt 



(A6) 
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From this dis cuss ion, it should be clear that th is fo rm is valid only for r oo. It is straightforward to see t hat in the 
same regime (A6) solves □$ = 0. Expression (A6) is the closest we can get to the usual plane wave form ([Al| ). The 



long-range gravitational field is simply taken into account by an appropriate phase modification. 

Next, we consider a plane wave in Kerr geometry. For simplicity we take the z-axis to coincide with the black hole's 
spin axis. We can then write 

^pianc « — - E c|''^^rV)^r (^)e'""^ (A7) 



The radial wavefunction is solution of 

, A + 2amuj f In 
+ + O 



dri 



= (A8) 



It is obvious that both (A7) and ( [A§| ) are different from the corresponding flat space expressions. That is, unlike in 
the Schwarzschild case, we are not able to derive an explicit form for a plane wave. Thus, we postulate the following 
asymptotic expression for a plane wave travelling along the z-axis 

^f. Aur-j, cos 9~iut / A ^^ 

•Pplanc = e (A9j 

where is the appropriate tortoise coordinate (|^). The field given by (|AS| ) is a solution of □$ = for r oo. For 
the general case of a plane wave travelling along a direction that makes an angle 7 with the z-axis the appropriate 
expression is given by (||). 

APPENDIX B: ASYMPTOTIC EXPANSION OF PLANE WAVES 

In this Appendix the asymptotic expansion of a plane wave in Kerr background is worked out. The calculation 
presented here is identical to the one found in the Appendix Al of jlj, but here it is specialised to s = 0. We have 
seen that for r — s- cx) the plane wave decomposition becomes 

g»^r.[sin7si„esi„^+co.7cose] _ J_ ^ ^(0) ^ (0) _ oo)^,"^ (6i)e""'^ (Bl) 

l,m 

where 7 is the angle between the wave's propagation direction and the positive z-axis. We multiply this expression 
by 5;°^, (6'')e^™ and integrate over the angles to get (after a trivial change I' ^ I, m' —f ni at the end) 

ChMllir 00) K Ujr d9 sine S{',^,{e)e''^'''''°''"=°'^ ^^^^<^r, sin-r sine sin ^-imv 

^0 "'0 
The integration over ip can be performed with a little help from , and the result is 

I.2TT 

dipe"^''* sin 7 sin 9 sin ^^i-m^ ^ ^^j^ (^^^ ^j^^ ^ gj^^ (gg) 

Since the Bessel function has a large argument it can be approximated as , 

V27rz ^ ' 
Note that this approximation is legal as long as 7 7^ 0. The o n-ax is case 7 = can be treated seperately, in a way 







similar to the one sketched here. Using this approximation in (B2) we get 



(B5) 

X±=/ de^/^SZ{^)e"^'''''°''^^^^^ (B6) 
Jo 

The T± integrals can be evaluated using the stationary phase approximation. We obtain 

<fy'lir ^ « 27r [(-^)™+le-'■-5^-(7) + z^+^e— ^f^l'T - 7)] (B7) 
We finally get ( |lO| ) by using the symmetry relation 

SZ{7:-6) = {-ir"^SZ{e) (B8) 
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APPENDIX C: CALCULATION OF SPHEROIDAL HARMONICS AND THEIR EIGENVALUES 



For the numerical calculation of the spheroidal harmonics we have adopted a "spectral decomposition" method, first 
developed by Hughes in the context of gravitational wave emission and radiation backreaction on particles orbiting 
rotating black holes. In the present work we have specialised this technique for the spin-0 spheroidal harmonics. The 
angular equation satisfied by S^^^{9) is, 



1 



de 



+ [{atoy cos^ e - 



sm 



Ebn]Si^ - 



(CI) 



where Eim denotes the corresponding eigenvalue. For the special case aw = we have Ei„i = l{l + l) and the solution 
of (|C1|) is the familiar spherical harmonic (here we are suppressing the dependence on ip) 



2l + l{l-m)\ 
Att {l + my. 



1/2 



Plni (cos I 



where Pim is the associated Legendre polynomial. It's numerical calculation is quite straigthforward 
the recurrence relation 



-Pim (a;) 



1 



[x{2l - l)Pl-l.rn -{l + m- l)Pl-2, 



I — m 

with "initial conditions" 

Pmmix) - (-l)"(2m - 1)!!(1 - 

P7n+l,rnix) = (2m + l)xPmrnix) 

We can always expand the spheroidal harmonic in terms of spherical harmonics, 



(C2) 
based on 

(C3) 



(C4) 
(C5) 



j=\m\ 



(C6) 



Substituting this spectral decomposition in (|Cl|), multiplying with Yim{0) and integrating over 9 we get 



where 



j=\m\ 



'o 



This integral can be evaluated in terms of Clebsch-Gordan coefficients 



2/ + 1 



< j2m0\lm >< j200|W > 



It follows that c"j ^ only for j -1,1,1 + 1. Then (C7) gives 



(C7) 



(C8) 



(C9) 



(CIO) 



We can rewrite ( plO ) as an eigenvalue problem for the matrix — {auj)^d^ with eigenvector 6* = bf^ and eigenvalue 
Elm ■ Clearly, M is a real band-diagonal matrix. Standard routines from can be employed to find the eigenvectors 
and eigenvalues of such a matrix. Then, the spheroidal harmonic is directly obtained from (CS) (even though it 
involves an infinite sum, in reality only few coefficients are significant). The described spectral decomposition 
method is reliable, unless auj becomes large compared to unity (under such conditions the matrix M is no longer 
diagonally dominant, and the convergence of the method is very slow). In effect, very high frequency cross sections 
for Kerr scattering will be inaccurate (especially when the black hole is rapidly rotating). 
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